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Abstract 

We prove the nonexistence of local self-similar solutions of the three dimensional in- 
compressible Navier-Stokes equations. The local self-similar solutions we consider here 
are different from the global self-similar solutions. The self-similar scaling is only valid 
in an inner core region which shrinks to a point dynamically as the time, t, approaches 
the singularity time, T. The solution outside the inner core region is assumed to be 
regular. Under the assumption that the local self-similar velocity profile converges to a 
limiting profile as t — ► T in L p for some p 6 (3, do), we prove that such local self-similar 
blow-up is not possible for any finite time. 



1 Introduction. 



In this paper, we study the 3D incompressible Navier-Stokes equations with unit vis- 
cosity and zero external force: 



ut + (u • V)u - 
V • u = 0, 
u\ t =o = u (x). 



-Vp + An, 



fl.l) 



We assume that the initial condition uq is divergence free and uq G L 2 
for some p £ (3, do). 



n 
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Many physicists and mathematicians have made a great deal of effort in under- 
standing the physical as well as the mathematical properties of the 3D incompressible 
Navier-Stokes equations. One of the long standing open questions is whether the so- 
lution of the 3D Navier-Stokes equations can develop a finite time singularity from 
a smooth initial condition PJ. Global existence and regularity of the Navier-Stokes 
equations have been known in two space dimensions for a long time 0. One of the 
main difficulties in obtaining the global regularity of the 3D Navier-Stokes equations is 
mainly due to the presence of the vortex stretching, which is absent for the 2D prob- 
lem. Under suitable smallness assumption on the initial condition, the local-in-time 
existence and regularity results have been obtained for some time [ZlElini. But these 
results do not give any hint on the question of global existence and regularity for the 
3D Navier-Stokes. 

In this paper, we prove the nonexistence of local self-similar singular solutions of 
the 3D Navier-Stokes equations. The local self-similar solutions we consider are very 
different from the global self-similar solutions considered by Leray (HJ. The self-similar 
scaling is only valid in an inner core region which shrinks to a point dynamically as the 
time, t, approaches the singularity time, T. The solution outside the inner core region 
is assumed to be regular and does not satisfy self-similar scaling. This type of local 
self-similar solution is developed dynamically, and has been observed in some numerical 
studies. Under the assumption that the local self-similar velocity profile converges to a 
limiting profile as t — ► T in LP for some p G (3, oo), we prove that such local self-similar 
blow-up is not possible. We remark that the nonexistence of global self-similar solutions 
has been proved by Necas, Ruzicka and Sverak in [HI]. The result of was further 
improved by Tsai in [TK] , 

We prove our main result by using a "Dynamic Singularity Rescaling" technique. 
This technique is simple but effective. Below we give a brief description of this tech- 
nique. Assume that the solution of the 3D Navier-Stokes develops a local self-similar 
singularity at x = at time T for the first time. A typical local self-similar singular 
solution is of the form 

u (*>*) = ^f= c/ (y> t )> pOM) = ^— t P(y,t), y = — p=, (1.2) 

for ^ t < T. We assume that u is smooth outside an inner core region {x, \x\ ^ 
(T — t) a } for some a > small. In particular, u(x, t) and p(x, t) are bounded for any 
fixed \x\ > as t — ► T. Using this condition, we can easily show that 

\U(y,t)\^C(T)/\y\, \P(y,t)\ < C{T)/\y\ 2 , for \y\ » 1, t ^ T. (1.3) 
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Thus, it is reasonable to assume that U S LP for some p £ (3, oo). But the LP norm of 
U may be unbounded for < p ^ 3. 

We assume that there exists a limiting profile U(y) £ LP as t — > T 

lim||^(t)-C7|| L p=0, (1.4) 

for some p satisfying 3 < p < oo. 

Next, we introduce the following rescaling in time: 

for ^ t < T. Note that by this time rescaling, we have transformed the original 
Navier-Stokes equations from [0, T) in t to [0, oo) in the new time variable r. Since u 
is smooth for < t < T, U is smooth for < r < oo. It is easy to derive the equivalent 
evolution equations for the rescaled velocity: 

U T + U + {y V)U + 2(17 • V)U = -2VP + 2AU, (1.6) 

with initial condition J7| r =o = yTuo^y / VT) , where U satisfies V • U = for all times. 
The problem on the possible local self-similar blowup of the Navier-Stokes equations 
is now converted to the problem on the large time behavior of the rescaled equations 
H1.6|) . By assumption (|1 .4(1 . we know that U — ► U as r — > oo in LP. We will prove that 
the limiting velocity profile actually satisfies the steady state equation of (jl.6p : 

U+ (y ■ V)U + 2(U- V)77 = -2VP + 2AC7, (1.7) 

for some P. Now it follows from the result of [H] that U = 0, which implies that 
]im T -> 00 \\U(t)\\lp = for some p £ (3, oo). 

The fact that lim^oo ||?7(r)||iP = is significant because it shows that the rescaled 
velocity field becomes small dynamically as r — > oo. It is easy to show that if the 
the solution U is small in the LP norm at some time, r m , the solution must decay 
exponentially in r for r ^ r m . The exponential decay in U gives a sharp dynamic 
growth estimate in terms of the original velocity field. In fact, it exactly cancels the 
dynamic singular rescaling factor, (y/T — t) , in the front of U. This gives us a uniform 
bound on the growth of LP for < t < T with p £ (3, oo), and consequently it rules 
out the possibility of a finite time blowup at T [TT | [T^ | Ifi], 

The nonexistence of local self-similar blowup of the 3D Navier-Stokes equations has 
some interesting implication. First, the assumption on the existence of a limiting self- 
similar profile, U, can be verified numerically if a local self-similar blow is observed 



3 



in a computation. Secondly, this result is related to a recent existence result by one 
of the authors [I] for the axisymmetric 3D Navier-Stokes equations with swirl. Let v r 
denote the radial component of the velocity field and r = \J x 2 + y 2 . The result in 
shows that if lim r _*o \ rv r \ = holds uniformly for ^ t ^ T, then the solution of 
the Navier-Stokes equations is regular for t ^ T. By the well-known Caffarelli-Kohn- 
Nirenberg result ^P, if the axisymmetric 3D Navier-Stokes equations develop a finite 
time singularity, the singularity must lie in the z axis. One of the most likely scenarios 
that would violate the condition, lim r ^o \rv r \ = 0, is the local self-similar blowup of 
the Navier-Stokes equations. The result presented in this paper would rule out such a 
possibility. For more discussions regarding other aspects of the Navier-Stokes equations, 
we refer the reader to [71 l2l IT31 Ifl] . 

The rest of the paper is organized as follows. In Section 2, we state our main 
theorem and present its proof. The proof is divided into four subsections. A couple of 
technical results are deferred to the appendices. 

2 The main result and regularity analysis. 

Theorem 1. Let u G L 2 (R 3 ) n L P (M. 3 ) for some p & (3,oo) and T be the first local 
self-similar singularity time. Assume that U(y,t) defined by converges to U in LP 
as t — > T . Then we must have T = +oo, i.e. there is no finite time local self-similar 
blowup for the 3D Navier-Stokes equations. 

Before we prove our main theorem, we state the following well-known (L^,LP)- 
estimates for the heat kernel in M 3 , e _ * A , where A is the Laplacian operator. 

\\e~ tA w\\ L t ^ c t~(%~V /2 \\w\\ L f, (2.1) 

||Ve _£A txj||i« < c r( 1+ ^H) /2 ||u>|| L p, (2.2) 

for 1 < p ^ q < oo, co depends on p and q only. In our analysis, we take q = p and 
p = p/2. For this particular choice of p and q, we can choose a constant, Co, such that 
the above two inequalities hold. Throughout the paper, we will use Co and c\ to denote 
various constants that do not depend on the individual functions, and use Cj (j = 1, 2) 
to denote various constants that depend on the initial condition, uq. We also define 

7 = 3/p. (2.3) 
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Since 3 < p < oo, we have < 7 < 1. 



Proof of Theorem 1. 

We will prove the theorem by contradiction. Suppose that T < +00. This means 
that the solution u to problem 111. Ill develops a singularity at t = T for the first time, 
but u is the unique smooth solution of 111. Ill for < t < T and is bounded in LP . 

We will divide the proof into four steps, which are given in the following four sub- 
sections. 



2.1 Dynamic singularity rescaling and a priori estimates. 

We make the following dynamic singularity rescaling of the 3D Navier-Stoke equations: 



1 T 

- log , 

2 & T-t' 



u(x, t) 
p{x,t) 



1 



Vr^t 



fU(y,r) 



(2.4) 



■P(y,r) 



for < t < T. 



Note that with this dynamic singularity rescaling, we transform the time interval from 
[0, T) in the original time variable t to [0, 00) in the rescaled time variable r. It is easy 
to derive an evolution equation for the rescale velocity field: 



' u T + U + (y • V)U + 2(U • V)U 

v-u = o, 

K U\ T=0 = VTuq (x) . 



-2VP + 2AC7, 



(2.5) 



Note that since u(x, t) is the unique smooth solution of the Navier-Stokes equations 
111. Ill for < t < T, U (x, r) is the unique smooth solution of the rescaled Navier-Stokes 
equations l|2.5p for < r < 00. 

Let 4>{y) = (<f>i,<fe,(f>3) be a smooth, compactly supported, divergence free vector 
field in M 3 and tP(t) be a smooth, compactly supported test function in (0, 1) satisfying 
Jq 1 ip(T)dr = 1. Multiplying (|2.5|) by ip(r — n)(j){y) and integrating over M 3 x [n,n + 1] 
for some n > 0, we obtain after integration by parts 



71+1 



(-ip T (f> ■ U + • U - VV • (0 (8) y) ■ U - 2^V(/> • (U <S> U)) dydr 



71+1 



i/jA(f) ■ Udydr, (2.6) 



5 



where ip is evaluated at r — n. 
By assumption ljl.4|) . we have 



lim ||i7(r) - U\\ LP = 0, (2.7) 

for some p > 3. Thus ||C/(r)||iP is bounded for r sufficiently large, and ||J7|U P is also 
bounded. Let U(t) = U + R(t). By l|2~7j) . we have lim^oo ||i?(r)|| LP = 0. Substituting 
U(t) = U + R(t) into (|2.6|) and letting n — > oo, we will show that all the terms involving 
R will go to zero as n —* oo. It is sufficient to prove this for the nonlinear term: 

n+l 

• (R R)dydr. 



ri 



Let q = p/(p — 2) > 1. Then we have 2/p + 1/g = 1. Using the Holder inequality, we 
obtain 



rn+l p r 

1/ / ^V<p- {R® R)dydr\ ^ C sup / \V(f)\\R\ 2 dy 

< C\\V(t>\\ Lq sup \\R(t)\\ 2 Lp ^ 0, as moo. 

Other terms can be proved similarly. Therefore, by letting n — > oo, we get 



\ T (r)dT) [ <P(y)U(y)dy 
o / Jm? 

+ (J ^(r)dr) (J J<f> -U - V • (0 ® y) -U - 2V0 • (U dy 



2 / ^(r)dr / A<j> ■ Udy . (2.8) 



o 

Since ip has compact support in [0, 1] , we conclude that 

/ ipr{r)dr = 0. 
Jo 

Moreover, we have ip(r)dT = 1 by assumption on ip. Thus, we obtain 



/ [<p ■ U - V • (<j> <g> y) ■ U - 2V<p -(U®U)- 2A<j) ■ U) dy = 0. 
it 3 



(2.9) 



Thus, U is a weak solution of the steady state rescaled Navier-Stokes equations: 

U + (y ■ V)U + 2(17 • V)7J = -2VP + 2AU, (2.10) 

with V • U = 0. Let Rj be a Riesz operator with Fourier symbol One can easily 

modify the proof of Lemma 3.1 of [10J to show that P = RjRk(UjUk)- 

Since U G L p for some p G (3, oo), we can apply Theorem 1 of [15J to conclude that 
U = 0. As a result, we obtain the following a priori decay estimate for ||E/"(t)||lp. 

Lemma 1. The solution U(x,t) of the rescaled Navier-Stokes equations 12. .5)) satisfies 
the following decay estimate: 

lim \\U(t)\\ lp =0. (2.11) 

T— >00 

For the purpose of our later analysis, we will choose a r m large enough to satisfy 
the following inequality: 

2c 2 oCl \\U(T m )\\ LP < i (2.12) 
where the constant ci is defined by 

9 1 1 + 7 

d = (l^ + 2)( 1 - e "T" a_ - ( 2 - 13 ) 
The reason for such a choice of r m will become clear later in our analysis. 



2.2 Dynamic decay estimates for the rescaled equations. 

In this subsection, we will perform estimates for the rescaled Navier-Stokes equations 
starting from r = r m with the initial value, U(x,T m ): 

' V T + V + (y • V)V + 2(V • V)V = -2VP + 2AV 
< V-V = (2.14) 

. V| r =o(a;) = Vo(x) = U(x,T m ), 

where r m is defined by (|2.12p - l|2.13p . Since U(x,r) is the unique smooth solution of the 
rescaled Navier-Stokes equations (|2,5|) for < r < oo, the function V(x,t) defined by 

V(x,t) = U(x,T + r m ), forr^O, (2.15) 

is the unique smooth solution of 1|2.14|) . 
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Next, we perform estimates for the linearized operator in 1)2.141 

dV 



+ V + (y ■ V y )V - 2A y V = 0, (2.16) 



with initial value V| T =o = Vq. 

Let y = e T y and V{y,r) = V(y,r). Then we have 



dV 

— + V - 2e~ 2r AyV = 0, (2.17) 
o- 



T 



with initial value V"| r= o = Vq. 

Taking the Fourier transform of 1)2. 17)1 . we get 

dV ~ ~ 

— + V + 2e- 2T \£\ 2 V = 0, (2.18) 

where the Fourier transformation is defined as /(£) = J f{x)e~ 2mx '^dx. Equation 1)2.18)1 
can be written as 

JL( e T + 2\eKe-*°ds^ {T) \ =0 (219) 



9r 

Integrating from to r, we get 

V(t) = e - T -l«l a ( 1 - e " ar )^. (2.20) 

Using the explicit formula of the Fourier transform of a Gaussian in three space dimen- 
sions (see, e.g. [TH]) 

e^N 2 = J-e-^l 2 /^ (2.21) 
cr 



with a 2 = f 1 _J-a T J , we obtain 

jr-l ( e -|CI 2 (i-e--)) = (—^\ 1 e-^M'/H-.*), (2.22) 

where J r " 1 /(^) = / f{i)e 2mx< di is the inverse Fourier transformation. Therefore, we 
have 

3 

V(y, t) = e~T (jZ^) 2 / V o (5) ( e — 2 l5--l 2 /(l-e- 2 -)) (2 . 23) 
Denote by e~ rA the solution operator of the linearized equations 1)2.17)1 . Define 

*o(t) = (1 - e~ 2r ), (2.24) 



and denote A as the Laplacian operator, then we have 

e — V = V(y, T ) = e~ T (e~ to ^ A V ) . (2.25) 

Define the following bilinear operator: 

F(U, V) = 2 (l - V (-A)" 1 V-) V • (U <g> V) . (2.26) 

In particular, if we set V = U, we have 

F(U,U) = 2(V ■ (U ®U) -V (-A) _1 V -V • (U ®U) 
= 2 (U • VU + VP) . 



(2.27) 



The rescaled 3D Navier-Stokes equations (|2.14)l can be converted into the following 
integral equation: 

V{t) = e~ rA V - ( T e-( T ~^ A F(U, U)(s)ds. (2.28) 
Jo 

To solve the integral equation ([2.28)1 . we construct a successive approximation, V^ n \ 
using the following iterative scheme (see [5]): =e- rA V , 

y(n+i) = y(0) _ G (y(n) iV W), n ^ 0, (2.29) 
where the bilinear operator G(U, V) is defined as follows: 

G(U,V)= [ T e- {T - s)A F{U,V)(s)ds. (2.30) 



To establish the convergence of the approximate solution sequence, V^ n \ we need 
to use the following lemma, which follows from 1)2.25)1 and the well-known (L q ,L p )- 
estimates (I2.1I)-(I2.2I) for the heat kernel. 



Lemma 2. Let V 6 LP for 1 < p ^ q < oo. We have 

\\e- TA V\\ Li ^ c e"( 1 - 3 ^>t (r)~(^"i) /2 ||y|| i p, (2.31) 



Ve- rA y|| Lq - < c e-( 2 - 3 / g >t (r) ^ /2 \\V\\ LP . (2.32) 



9 



The lemma can be proved easily by noting that the heat kernel actually acts on the 
variable y through the function V(y, r) and y = e~ T y. Thus we lose a factor e 3r ^ when 
we estimate the L q norm by changing variables from y to y, but we gain a factor of e~ T 
when we differentiate with respect to y. 

Applying ()2 . 3 1 f) with p = q = p, we obtain 

\\V^\\ LP (r) = \\e~ rA V Q \\ L v(T) < coe-^lVoIlL*, (2.33) 

where 7 = 3/p. To estimate V)[|£p, we use (|2.H2|I with q = p and p = p/2: 

1+7 



\\G(U,V)\\ lp (t)^2co / e-( 2 -^^- s )t (r- S )-^||C/|| L p( S )||F|| L p( S )d S , (2.34) 
J 

where we have used the Holder inequality \\U (8> V^|| iP /2 ^ H^IIlp H^Hl? and the fact 
that (-A)" 1 V-V- is a Rietz operator of degree zero, which is a bounded operator from 
LP to LP. In particular, we obtain by setting V = U that 

\\G(U,U)\\lp(t) ^2c [ T e-( 2 -^ T - s h (T - sy^imUs^s . (2.35) 
Jo 

Now, applying l|2.ririj) and (|2.rif)|l to the iterative scheme l|2.29j) . we get 

||^ {n+1) ||L P (r) s$ c e-^ T \\Vo\\ L p + 2c f e'^^to (r - s)^ \\V^\\ 2 LP (s)ds. 

Jo 

(2.36) 

Define 

K n = sup \\e^ )T V^(r)\\ LP . (2.37) 

0^T<OO 

Multiplying l|2.rif)|) by e^^ T on both sides and using l|2.37|) . we obtain 

e^-^WV^W L P (t) < co\\Vo\\ LP + 2c e- T KlJ\^ (l - e" 2 ^) ds. (2.38) 
In Appendix II, we will prove that 



7»h _ p- 2 ( r - s ) 



1+7 

1 _ e -^ r -*M 2 j s ^ Cl; forallr^O, (2.39) 

where ci is defined in l|2.1rij) . Now, take the supremum of the both sides of I|2.ri8j) for 
all r ^ 0, we obtain the following recurrence inequalities: 

K n+l s$ K + MR 2 , for n > 0, (2.40) 

with K n \ n= Q = Kq : where 

M = 2c ci, K = c \\Vo\\lp- (2.41) 
We will prove the following lemma in Appendix I. 
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Lemma 3. Let Kq and M be two positive constants satisfying 

K M < i, (2.42) 

D 

then there exists a positive constant K max , such that 

K n < K max , for all n ^ 1 (2.43) 

ZioWs /or i/ie recurrence sequence K n satisfying \2.4Uj) - Moreover the upper bound K ma/X 
satisfies 

2MK m , x ^ 1 (2.44) 
Recall that in (|2.12|) . we have chosen r m such that 

2cgci||Z7(r m )|| L p <l/6. (2.45) 

Therefore, we have 

K M < 2cg ci ||C/(r m )|k P < J- (2.46) 

o 

Thus, for our choice of r m defined in l|2.12jl . the recurrence sequence K n has an upper 
bound K max for all n. That is 

\\V^\\ LP (r)^K max e-^>, for n > 1. (2.47) 

2.3 Convergence of the approximate solution sequence. 

In this subsection, we will establish the convergence of the approximate solution se- 
quence, {T^ n )}, and study the property of its limiting solution. We will first show that 
the approximate solution sequence {y( n )} is a Cauchy sequence in LP . By subtracting 
l|2.29|) with index n from that with index n — 1, we obtain 

\\V {n+1) -V (n) \\ LP 

= \\G(V( n \v^) - G(Vl n - 1 \v( n -V)\\ Lt , 

= \\G(V {n) ,V {n) - V {n ~ l) ) + G(V {n) - V {n - l \ V^ n - X) )\\ L P. 
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Using (E3H), ll2~37ll and fftTSfy . we obtain 

e (l-7)T||^(n+l) -V^W^t) 

< 2c e (1 - 7)T J\- {2 -^ T - s k (r - s)^ (||U (n) || LP + H^" 1 ) || LP ) - T^" 1 ) 

^4c K max e- T [ T e^toir-sy^dsf sup e^'H^W - \\ LP (s)) 

JO \0^s<oo / 

^4c 0Cl ^ max sup e^ s \\V^-V^\\ LP (s) 

0^s<oo 

sup e (1 - 7)s ||y (n) -U^-^Ulp^), 

2 0^s<oo 

where we have used H2.41|) and l|2.44fl in deriving the last inequality. Taking the supre- 
mum on the left hand side would yield 



sup ,, ... 

0^T<OO ^ 0^T<OO 



3 (l-7)T||y(n+l) _yW|| Lp ( r ) <; * gup e (l-7)r|| F (n)_ y (r l -l)|| Lp ( r ) j (2>48) 
which implies 

sup e ( 1 -^||y("+ m ) - y (n) || LP (r) < d (\ ] , for any n,m > 1, (2.49) 



0<T 



where Ci depends on only. Thus {lA n )} is a Cauchy sequence in BC ([0, oo) ; LP (K 3 )) 
Here BC ([0, oo) ; L p (M 3 )) denotes the class of bounded and continuous function from 
[0, oo) to L P (M. 3 ). As a result, we have proved that (r) converges uniformly to a 
limiting function V (r) in BC ([0, oo) ; L p (R 3 )) . Taking the limit n -> oo in (|2~I7|) . we 
obtain 

||V|| LP (r) ^ ^ max e- (1 - 7)r . (2.50) 



Next, we will show that V is a solution of the integral equation l|2.28p . To this end, 
we define R^\x,t) = (x,t) — V(x,t). We have just shown that 

sup e (1 -~ t)T \\R( n) \\ LP {T) = sup e {1 -~ t)T \\V {n) -V\\ L p{t) -» 0, (2.51) 

0^r<oo 0^r<oo 

as n — > oo. Now substituting = V + i?( n ) into the iterative scheme l|2.29|) and 
using the bilinearity of operator G(U,V), we get 

V - V {0) + G(V,V) = -{R {n+1) + G{R {n \V) + G(V, i? (n) ) + G{R {n \ R {n) )). (2.52) 
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We will prove that the error terms on the right hand side of (|2.52|l tend to zero uniformly 
for all r ^ 0. It is obvious that ||i?( n+1 ) \\lp — * uniformly as n — » oo from 1)2.51)1 . 

To show that the error terms which are linear in RS n ) tend to zero uniformly, we 
use 1)2.34)1 and the a priori bound on V given by 1)2.50)1 . Specifically, we have 

\\G(R^,V) + G(V,R^)\\ lp (t) 

^4c [ T e-( 2 -y^-s) to ( T _ s y±P \\R(n)\\ LP ( s )\\V\\ LP (s)ds 
Jo 

^ 4c ^ max e-( 1 -^e- T f eiH {r-s)-^ds( sup e^ s \\R^\\ LP (s)) 

JO \0^s<oo / 

< 4c c 1 e-( 1 -^^ max sup e^ s \\R^\\ LP (s) 

0^s<oo 

< sup e ( - 1 " r >\\R^\\ LP (s)^0, 

0sCs<oo 

uniformly for all r as n — ► oo, where we have used M = 2cqC\ and 1)2.44)1 . 

To show that the nonlinear error term G{RS n \R> n ') also tends to zero uniformly, 
we note that the a priori bounds on and V also provide the following a priori 

bound for R^: 

\\R^\\ LP (r)^2K max e-^ T , for n > 1. (2.53) 
Using 1)2.53)1 and applying the same argument as above, we can prove that 
\\G{R^\rW)\\ lp {t) ^ sup e ( 1 -7)«||i J (»)|| £l ,_o j 

0^s<oo 

uniformly for ^ r < oo as n — * oo. 

Now, passing the limit n — ► oo in the LP norm, we obtain 

V(t) = V {0) - G(V,V), for all r > 0, (2.54) 

which shows that V is a solution of the integral equation l)2.28p , satisfying the decay 
property 1)2.50)1 . 

2.4 The non-blowup estimates in the original variables. 

In this subsection, we will complete the regularity analysis in the original physical 
variable. By the uniqueness of strong solutions in L p with p > 3, we have 

\\V\\ LP (r) = \\U\\ L p(T + T m ), for0sSr<oo. (2.55) 
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Now we can use the decay estimate for V in Ij2.50p to obtain a decay estimate for 
U, which in turn will rule out the possibility of a finite time singularity for the 3D 
Navier- Stokes equations. 

Using (|2,5()p and 1)2. 55 j) . we immediately obtain a decay estimate for U: 

\\U\\lp (t) < WNK), for r > r m . (2.56) 
This proves the following decay estimate for U. 

Lemma 4. The solution U(x,t) of reseated Navier-Stokes equations 12. .5)) with T m 
defined by \2.12\) has a uniform decay rate in r as follows: 

\\U\\ L v (t) ^ K max e'^-^ T - T ^ , for r > r m . (2.57) 

Substituting the relation 

u(x,t) = -^=U(y,r) (2.58) 

into (f2"57) . we obtain for t m sC t < T with t m = T (l - e _2T " m ), 
(T - tV/ 2 

\w\\Mt) = \ t _; )1/2 \\u\\l*{t) 

< K max (!_ 7 )( T _ Tm ) _ ^max^J^_ (i_ 7 ) r 

^ (r_f)(l-7)/2 (T-i)(l"7)/2 

x™,^ 1 -^ /t-A (1 ~ 7)/2 



(T-t)(i-7)/2 ^ T 

K p( 1 -7)T"m 

xl max c 

T(l-7)/2 



(2.59) 



<^ t^^—. , for t m ^i<T. 



Since uo G £ p for some p G (3, oo), it is easy to show that there is a local-in-time 
smooth solution whose L p norm is bounded [2] (This can also be proved directly by 
using the same iterative scheme applied to the original Navier-Stokes equations for a 
short time). Moreover, since T is the first singularity time, we conclude that u is smooth 
for < t ^ t m < T and has a bounded LP norm for t ^ t m . Thus, ||u||i>(i) is uniformly 
bounded for ^ t < T. 

Now, we can apply the so-called Ladyzhenskaya-Prodi-Serrin condition (see [6], 
and [12J), which is also known as the L p ' q criteria. The so-called L p,q criteria state that 
if a suitable weak solution of Ijl.ip satisfies 



u 



G L q ([0, T); LP (K 3 )) (2.60) 
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with ^ 

- + p€[3,oo], (2.61) 

p q 

then u is a smooth solution of the 3D Navier-Stokes equation up to t = T. In our case, 
we have obtained a uniform bound in IP for u with p £ (3, do) for ^ t < T. Thus 
the L p ' g criterion is satisfied with q = oo. Therefore, we conclude that u is a smooth 
function in M 3 x (0,T]. 

This conclusion contradicts with our assumption that u would cease to be regular 
at time T for the first time. This contradiction implies that u can not develop a local 
self-similar singularity in any finite time. This completes the proof of Theorem 1. 

Appendix I. 

In this appendix, we prove Lemma El 
Proof of Lemma El It is sufficient to obtain an upper bound for the recurrence 
equalities 

K n+1 = K + MKl K = K . (2.62) 

It is easy to see that K n ^ K n , for all n ^ 1. To simplify the notation, we will drop 
the tilde in K n in the following. Define l n = K n+ \ — K n , then we have 

l n = M (K n -i + K n ) k-i. (2.63) 

It is easy to see that l n > for all n ^ and K ri is a monotonely increasing sequence. 
We claim that 

M {Kj-i + Kj) ^ -, for all j > 1. (2.64) 

We will prove (|2.64|1 by an induction argument. 

1. For j = 1, we have 

MiKo + K^ = M(K + K + MK$) < i (2.65) 

from the assumption KqM ^ |. 

2. Assume that Ij2.64|l holds for all j ^ n, we will prove that it also hold for j = n+ 1. 
Let a = 1/2. It follows from l|2.63|l and the induction assumption that 

lj ^ alj-i, for all 1 ^ j ^ n, (2.66) 

which implies that 

lj < oPIq- (2-67) 
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Applying K n+ \ = K n + l n recursively and using (|2.67p . we obtain 



^ K + /oEi=o« J 
1 - a n+1 

= Kq + Iq— 

1 — a 



(2.68) 



where we have used MKq ^1/6. Define K mSLX = ^Kq. Then we have 

2MK max = ^MK ^l<\- (2-69) 

Thus, we obtain 

M (K n + K n+X ) 2MK max < 1 (2.70) 
This proves the claim 1)2.64(1 by induction, and we obtain 

K n ^ K mSLX , for all n > 0. (2.71) 

We have already shown that 2MK max < — in 1)2.69(1 . This completes the proof of 
Lemma 



Appendix II. Proof of estimate (12.391) 



In this appendix, we prove estimate 1)2.39(1 . First, we state a useful inequality 

|1 - e~ 2x \ > (1 - e~ 2 )|x|, fbrO<i<l, (2.72) 

which is a consequence of the fact that (1 — e~ x )/x is a monotonely decreasing function 
for x > 0. We consider two cases. If r > 1, we divide the integral into two parts as 
follows: 

. 1 + 7 

J r eT s (l-e- 2 ( T - s )) 2 ds 



= Ic ,+J^i« 7 '(l-«- 2(T - ) )" 2 ds 
< j;- 1 e7S (1 " e " 2 ) _i? * + /.-I « 7T ((1 - e" 2 )(r - S ))-^ * 

= (l-e- 2 )--^^-e-(l-e- 2 )-^ T ^ ( r- S ) i - 

1+2 ^<v-n -, . _ , n . _ 1+2 



S = T 



(2.73) 



T-l 



7 1-7 

_ Wi _ „- 

1-7 7, 
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where we have used ()2.72[1 . Thus we prove 




1 + 7 



e 7« (l _ e -2(r-s)\ 2 ^ ^ Cie -(l- 7 )r < 



for all r > 1, 



where c\ is defined in (|2 . 1 3|) - For r ^ 1, we have by using (|2.72p 




1+7 



/ 

Jo 



T 



((1 - e~ 2 ) (r - s)) ^ ds 



— e 



2(t-«) 



ds 



e 



(l- 7 )r 



2 e -( 1 -T) T 




1-7 



1-7 



( X -^ T ^ ci. 



(2.74) 



This proves (l2~^?ll . 
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